Error characterization in Quantum Information Processing: a protocol for analyzing 
spatial correlations and its experimental implementation 



Cecilia C. Lopez, Benjamin Levi, David G. Cory 
Department of Nuclear Science and Engineering, MIT, Cambridge, MA 02139, USA 

We present a protocol for error characterization and its experimental implementation with 4 qubits 
in liquid state NMR. The method is designed to retrieve information about spatial correlations and 
scales as 0(n w ), where w is the maximum number of qubits that have non-negligible interaction. 
We discuss the practical aspects regarding accuracy and implementation. 
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PACS numbers: 03.65.Yz, 03.67.Ac, 33.25.+k 



I. INTRODUCTION 



Precise and reliable control of the system constituting 
a quantum information processor (QIP) remains one of 
the biggest challenges in the quantum information field. 
In order to assess the reliability of a device and to tai- 
lor quantum error correction schemes to a faulty one, 
we need to characterize the errors occurring in the sys- 
tem. Quantum Process Tomography (QPT) [l[ presented 
a first answer to this problem, providing full character- 
ization of the process under analysis. Experimental im- 
plementations of QPT have been already conducted in a 
variety of small systems [H, El, 0, [5j] . Nevertheless, QPT 
becomes impractical beyond a few qubits, as it requires 
0(2 4n ) experiments for a system of n qubits. In recent 
years, the idea of getting less information at a lower cost 
has become a popular strategy to tackle error charac- 
terization, and several works have been devoted to the 
subject 0, 0, @, ©, [HJ [0, [3. Our proposal fits in this 
context, providing a subset of information yet still using 
scalable resources. 

The scheme we present selectively keeps information 
about the spatial correlations of the errors occurring in 
the process under study (a gate, a noisy channel, etc.). 
Both the magnitude and the structure of the errors are 
relevant to evaluate fault-tolerance. In particular, fault- 
tolerance threshold theorems are designed for certain 
conditions of spatial correlation (also termed range or 
locality) [HI]. So even once it is experimentally deter- 
mined that only up to w qubits are involved in an error 
process, we have to further establish in which way the 
(") possible sets are being affected. 

We have implemented our protocol in a liquid state 
NMR 4-qubit QIP. The core mathematical work for this 
protocol was introduced in Q ■ Here we extend our pro- 
posal to a more general setting and include an experi- 
mental realization. Our basic method, like others pro- 
posed HUSH 

[ToL [Til ] , assumed error-free implementa- 
tion stages. This idea is of course unrealistic in practice, 
and implementation errors complicate the task of reli- 
able error characterization. Thus here we have included 
an analysis of their effect. 



II. THEORY BEHIND THE PROTOCOL 

We start by describing the action of a general map on 
the state of an n-qubit system (D — 2"), described by 
an initial state po, as 



S(po) = / P(v) E(ff) p E\ff) dff. 



(1) 



The vector fj denotes D 2 complex coefficients {770, r)i, I = 
1, . . . , D 2 — 1} that parametrize E, an arbitrary operator 
in the Hilbert space Hd, as 



E = r l0 I+ VI O u Oi =(g)0 



(./) 



(2) 



1=1 



where each O] is an element of the Pauli group 
{I,a Xl a y ,a z }, but at least one factor in each Oi is a 
Pauli matrix. Notice that Tr[0;0/'] = D5ij>. When p(fj) 
is a (real) nonnegative distribution, eq. |1| describes an 
arbitrary completely positive (CP) map, and the con- 
dition J p(rf)\fj\ 2 dfj — 1 guarantees the preservation of 
Tr[p]. This representation of a CP map is an operator 
sum representation with continuous parameters. We pre- 
fer this form as it's more suitable to describe non- unitary 
dynamics arising from stochastic Hamiltonians. Our pa- 
rameters, the 77, are trivially related to the parametriza- 
tions of the operator sum representations used in other 
works 

[1G3. Also, for small r/i, it's possible to relate ([T]) 
to a description of the noise in terms of generators rather 
than propagators 0]. 

When necessary, we shall denote the 77/ in more 
detail as rfj'k", where j > k > ... label qubits, and 
p,q,... = x,y,z. Therefore rf^" labels a term in eq. 
([2]) that is a product of <j p for qubit j, a q for qubit 
k, etc., and / for the qubits absent in the subscript. 
Notice that the number of qubits in the subscript gives 
the Pauli weight (also Hamming weight) of the term [141 ]. 

In our protocol we measure a subset of m qubits at 
a time, which will allow us to extract the magnitude of 
the errors involving that subset. So we now break the 
system in two: m qubits belonging to the Hilbert space 
Hm, which will be the ones measured, and m = n — m 
qubits belonging to the complementary space Hjt, H-d = 



2 



Hm ® "Hp-. We require the initial state to be separable 
in these two spaces, and within TLm, thus 



Po — P ® Po 



and Pq M ^ 



(3) 



We now perform a £/(2)® m twirl on the target map S, 



(4) 



k=l 



where the Ck are m-fold tensor products of the twirl op- 
erators on one qubit, for the m qubits in Hm- For the 
task we propose, a subset of K = 6 operators from the 
Clifford group will suffice. We refer to this minimum set 
of operators required to apply the twirl (0| as the 6 m - 
Clifford element pool (see Appendix A for references and 
some mathematical details on this twirling). 

The reduced density matrix p^ = Tr^j[pi] is the 
state of the m qubits we measure. In general, the fi- 
delity decay, which gives a measure of how p\ differs 
from the initial po, encodes information about the map 
S we are characterizing [y, |7|. In particular, consid- 
ering only the qubits we will measure and expressing 



(M) 

Po 



(j®m + p ^ v >)/2 m , we obtain 
Tr[U M) ) 2 ] - ^m[ P [ M) P [ M) ] = 7 (M) - n&> (5) 



v( M ) 



(6) 



jeM 



and fiW = if p%> = for p^> ± can be 

exactly computed if desired). The quantity defined in 
the LHS of (SJ is the fidelity decay rate 0, 0] we will 
analyze in this work. In eq. © we have denoted (...) = 

J p(rf) . . . drf, Vj = Tt[(pq^) 2 ] the initial purity of each of 
the M-qubits, and 



,(A?) 



)(M) 



(2/3)(l-V j /2) if Of =a x ,a y ,a z 
V, if 0\ j) = I 



(7) 



For the derivation of eqs. ((MSI) we use the equivalence 
between a Clifford twirl and a Haar twirl [lH ] , and apply 
the tools developed in See also Appendix A. 

In the sum ([S]), a Z-term vanishes when Oi is the iden- 
tity operator for the m qubits being measured. So by 
choosing different sets of M-qubits, it is possible to leave 
out certain rji in a given 7. Note however that Cj does not 
distinguish the direction of the Pauli matrices, so there 
is an implicit coarse-graining of all the Oi that have the 
identity / for the same subset of m qubits. This is how 
we get the "collective coefficients" rjf° l , 



I'/ 



oil 2 



\V 



,PI 2 



etc. (8) 



For example, the terms and contribute 

collectively to 77^3 . 

By combining the 7's from different M-sets it is pos- 
sible to further isolate the collective coefficients. If we 
prepare a and 6 in a pure state, so V a = Vh = 1, we 
obtain 



!( 7 c 



U 7 0»_ 7 (<M>)) = (|^|») +J2(\<U 2 ) + ■ (9) 



(Please refer to Appendix B for the formulae of the 7's in- 
volved in this example). If 3-body and higher multi-body 
terms can be neglected, eq. © gives (h^l 2 )- Similarly, 
the combination of the seven 7W, 7 ^ J > and 7^ J J ^ for 
a set of three qubits a,b,c would return (\r]^ l c | 2 ), and so 
on. The collective coefficients then report on the spatial 
correlations of errors. 



III. THE PROTOCOL 

We present now a systematic protocol for measuring 
the collective coefficients involving the m qubits of a par- 
ticular subset. 

i. Prepare each of the qubits to be measured in the 

initial state |0). Prepare each of the other m qubits 

in the maximally mixed state 1/2. 
ii. Apply one of the m-fold Clifford operators from the 

6 m -Clifford element pool. 
Hi. Implement the target gate or noise under study. 
iv. Invert the Clifford operator applied in ii. 
v. Measure the projection of the resulting state on 

the initial state |0), for each of the qubits being 

measured. 

To implement the twirl, we repeat i — v each time taking 
a different operator from the 6 m -Clifford element pool, 
and average the results. 

In a canonical QIP, the implementation of this protocol 
to measure the decay rates involving m qubits will require 
N realizations. This will take care of: 1) preparing the 
desired initial state as in step i (starting from the 10)®" 
state, and randomly flipping the m qubits we don't mea- 
sure); 2) measuring the fidelity decay through repeated 
projective measurements of the m qubits as prescripted in 

(AT) JM), 



step v (notice eq. ([5]) becomes I — Tr m[Pq Pi ] = 7 



p—x,y,z 



for the proposed initial state); 3) implementing the twirl 
approximately, by randomly drawing the twirl operators 
for steps ii and iv from the corresponding 6 m -Clifford el- 
ement pool (or from an infinite pool of one-qubit random 
rotations). 

With this strategy, the outcome of the measurement 
step v is a Bernoulli variable, and thus N can be esti- 
mated from usual statistics. More precisely, the stan- 
dard error in our estimation of the decay rate will be 
<7 7 < l/\/~N (following the Central Limit Theorem), so 
for a desired precision S we must have N > S~ 2 . On 



3 



the other hand, the Chernoff Bound tells us that for a 
desired precision S and an error probability epf ^ 1, we 
must have N = log(2/e,/v)/(2(5 2 ), which is a stronger re- 
quirement when 6 < 2e~ 2 . In any case, N is independent 
of the number of qubits thus the efficiency of the protocol 
is independent of the size of the system. 

In the case of liquid state NMR ensemble QIP, pseudo- 
pure state preparation allows for initialization in the 
1/2 state over the ensemble of molecules, and ensemble 
measurements avoid the need of repeated realizations in 
order to perform step v by Quantum State Tomography. 
This is the case in our experiment. 

The measurement step v retrieves the information to 
calculate the decay rates for the m qubits and for any 
smaller subset of qubits within those. For example, 
twirling qubits a and b only (m = 2), we can obtain 
ry( a \ -yW and 7( Q ' b ) in one shot, and calculate (|?7™h| 2 ) as 
in eq. ([9]) neglecting the higher multi-body terms. This 
procedure can be repeated for the (™) = n(n — l)/2 pairs 
of qubits, and by doing so all the collective coefficients 
for 1- and 2-body terms can be extracted. 

The scalability of the method goes as follows. If 
we can neglect the multi-body terms above a certain 
Pauli weight w, and N is the number of realizations 
required to measure the fidelity 7 for w qubits, then 
with N rM < Nn w /w\ experiments we can estimate 
all the non-negligible coefficients. This should be com- 
pared against N2 4n , the overhead in the number of ex- 
periments required for QPT. We emphasize here that our 
proposal seeks to characterize the correlations among up 
to w qubits in order to establish the range of the noise - 
not to characterize the process fully. 

To use our protocol, the negligibility of multi-body 
terms above a certain Pauli weight w must be estab- 
lished a priori. In a canonical QIP we could apply our 
protocol to measure all the n qubits, obtain all the decay 
rates 7^, 7^'^, . . . ^C 1 '---'™) and extract all the collec- 
tive coefficients after only N experiments (independently 
of n). In this way we can handle all the Pauli weights, 
from to n. But the error in the decay rates ct 7 propa- 
gates into the rf° l for m qubits inefficiently, roughly as 

a v 2 cx ^Jj2j=o (j)°V Therefore our strategy of looking 
for spatial correlations after establishing a cut-off Pauli 
weight w: to our knowledge, neither QPT nor other pro- 
posals so far are able to make use of the negligibility of 
high order correlations in order to gain further insight. 

An approach to establishing this cut-off w, demand- 
ing the same resources as our protocol, is to apply the 
method developed by Emerson et al. [8J which gives the 
probability of errors happening distinguishing them only 
by Pauli weight (this is, an average of all the (™) col- 
lective coefficients having Pauli weight w). It is worth 
pointing out that both methods require the same exper- 
imental work: the algorithms retrieve different informa- 
tion because the measurement and processing of the data 
is different, but actually both protocols can be used com- 



plementarily. 

IV. EXPERIMENTAL RESULTS 

We implemented our protocol in a liquid state NMR 
QIP using the four 13 C-labeled carbons of crotonic acid 
[17| in a 400MHz Bruker spectrometer. The initial state 
preparation and all the gates required by the protocol 
were implemented with RF pulse sequences engineered 
using either GRAPE [lj] or SMP [13(c)] methods. Their 
simulated gate fidelities F g were on average 0.98. The 
typical experimental performance of one-qubit gates on 
the spectrometer is 1-2% below their simulated fidelities. 
See [2j for details on the model used in the simulation. 

We studied the following processes (see Table U for 
more details): i) A time suspension sequence Ie, since 
it is important to study our ability to "do nothing" in a 
system with a natural Hamiltonian that is always on. ii) 
An engineered error creating a coupling between qubits 
1 and 2, of the form C^iP) = exp(— z/3<r^cr 2 ). We chose 
[3 = 0.1, and (3 = 0.4 (the previous one applied consecu- 
tively 4 times). Hi) A CNOT gate between qubits 1 and 

2: CNOT = 0.5(7 + <t[ 1} + - a^a^). Also, this 
gate applied twice: CNOT 2 = I. These gates are more 
complex than one-qubit operations (which are typically 
less than 1 msec long) and they all involve refocusing 
idle times (periods of free evolution under the internal 
Hamiltonian) in their pulse sequences. 

The results on the measurement of the collective 
coefficients for the qubit pairs (1,2), (2,3) and (1,4) are 
presented in Table |U The pair (1,2) is the one targeted 
by the C12 and CNOT gates, while the other two are the 
pairs involving qubits 1 or 2 with the highest J-coupling 
[Tn |. We expect the errors for the three chosen pairs to 
be larger (due to internal evolution that is not perfectly 
refocused). See Appendix C for more details on the 
experiment and simulations using liquid state NMR QIP. 

The resulting (|?7™b| 2 ) for the various pairs of qubits 
shown in Table |T] exhibit good agreement with the pre- 
dicted ones. Notice that the largest differences between 
measured and predicted appear on the pair (1,2) and 
on the most complex gates: CNOT 2 , CNOT, C 12 (0A). 
This indicates that these deviations are due to the errors 
expected from spurious processes in our QIP, particularly 
an imperfect refocusing during the gate sequence, rather 
than from an imperfect implementation of the protocol. 

The simulations account for well-known sources of er- 
ror in liquid NMR QIP (imperfect pulse design [13(c)], 
RF field inhomogeneities [2j], and the presence of the 
magnetically active hydrogens in Crotonic Acid [17]). It 
worth mentioning that the main errors ocurring in one- 
qubit gates are correlated one-body errors (this is, a one- 
qubit rotation that is slighly off), and they do not in- 
troduce two-body errors, which are the main target of a 
spatial correlation analysis. This can be noticed also in 
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Gate 




<l^1| 2 > 


<ICil 2 > 


Ie 

iz.z msec 


meas 


0.02 


0.02 


0.01 


theo 


0.00 


0.00 


0.00 


F g = 0.96 


simG 


0.01 


0.02 


0.00 




simE 


0.01 


0.03 


0.00 


Cl2(0.1) 

4.88 msec 


meas 
theo 


0.02 
0.01 


0.02 
0.00 


0.01 
0.00 


F g = 0.99 


simG 


0.02 


0.00 


0.00 




simE 


0.02 


0.00 


0.00 


Ci2(0.4) 
19.52 msec 


meas 
theo 


0.26 
0.15 


0.03 
0.00 


0.03 
0.00 


Fg — 0.87 


simG 


0.23 


0.01 


0.00 




simE 


0.24 


0.02 


0.02 


CNOT 
11.00 msec 


meas 
theo 


0.32 
0.25 


0.01 
0.00 


0.03 
0.00 


F 9 = 0.99 


simG 


0.25 


0.01 


0.00 




simE 


0.28 


0.02 


0.00 


CNOT 2 
23.76 msec 


meas 
theo 


0.07 
0.00 


0.05 
0.00 


0.04 
0.00 


F 9 = 0.97 


simG 


0.01 


0.01 


0.00 




simE 


0.01 


0.02 


0.00 



TABLE I: Measured (meas) collective coefficients for selected 
pairs of qubits, for the various gates studied, theo are the the- 
oretical values for these gates as described in the text. The 
sim values come from numerical simulation: simG are the 
values obtained by the calculation of the propagator as given 
from the RF pulse sequence corresponding to the gate (G) 
alone; simE are the values obtained by completely mimick- 
ing the experiment (E), considering the average of 36 simu- 
lations, each applying RF pulses to prepare the initial state 
and implement the Cliffords gates and the gate under study. 



the fact that the simG values are similar to the simE 
values (a change of 0.00 — 0.01, except in two cases where 
we found 0.02 and 0.03). Moreover, there is a contribu- 
tion arising from Ti relaxation. However, a calculation 
over the theoretical propagator and the numerics over the 
simulated one show a change on the order of 0.01. 

Elements outside our model of the system, which 
would explain futher the gap between theory and 
experiment, are Bq (static) field inhomogeneities, the 
presence of transients and residual non-linearities in the 
spectrometer circuitry, and an imperfect spectral fitting 
of the measured signal. These are well-know issues in 
liquid NMR QIP, whose effect falls within the 1 - 2% 
error. 

We must differentiate between the implementation er- 
rors in the protocol (initial state preparation, one-qubit 
twirling and readout), and the errors in the gate under 
study. The former ones affect the accuracy of proto- 
col, which accounts for the measured non-null coefficients 
that are expected to be zero. As discussed, RF field inho- 



mogeneities, the presence of hydrogens and T2 relaxation 
already give an error bar 0^2 w 0.03. And there are still 
other sources of error mentioned above that could make 
a v 2 larger, but within that order. 

Moreover, a fiducial initial state preparation is critical 
to the sucess of the algorithm. We can quantify this 
as follows. If we call £0 the error in the initial state 
preparation, and similarly we call E\ the error in the im- 
plementation of the Clifford gates, an error propagation 
in the formula for the 7's gives cr 2 < £q(1 + 47) + e\ , 
and then it's simple to propagate this into the formulas 
for the |t7 ; 2 |: for example, for a pair of qubits a, b we 

follow eq. (HI) and obtain oy = | J~tf a + tf b + cr^ & . 

These e's account for non-statistical errors (typically 
correlated one-body errors) and set the accuracy of the 
method. Given the low complexity of initial state prepa- 
ration and one-qubit operators, these e's are smaller 
than the errors in target operations (a fact reflected, 
for example, in the gate fidelities). This is why even 
though the theory was developed for error-free initial 
state preparation and twirl, the actual implementation 
can still retrieve information about the target operations. 



V. DISCUSSION ON THIS AND OTHER 
PROPOSALS 

Our protocol belongs to the family of characterization 
strategies based on the use of twirling to coarse-grain the 
original 0(2 4 ™) complex parameters to a poly(n) num- 
ber (c/ the work by Emerson et al. in [g, |8|). In our 
case, we gain detail about the process under study (spa- 
tial correlations). Note also that the required resources 
in our proposal are within the minimal performance ex- 
pected from a functional QIP: fiducial state preparation 
and readout, and a set of 6n one-qubit gates. 

One of the protocols presented by Bendersky et al. 
in [Tol | also returns similar information (although our 
method performs a coarse-graining of the directions). 
Nevertheless, their proposal requires more demanding re- 
sources (although still scaling as poly(n)): a full twirl on 
U(D) and the implementation of the Oi operators. 

Other characterization methods include the ancilla- 
assisted ones in [ToL HH , but unfortunately they have a 
rather strong requirement: one or more clean error- free 
qubits within the system. Contrasting with all the 
proposals discussed so far, the method developed by 
Knill et al. in [9] does not require error-free stages, 
allowing for certain type of errors to occur during 
the whole computation. Unfortunately it's not yet evi- 
dent how to take their scheme beyond one-qubit QIP [12] • 

On a different note, we would like to point out a 
particular feature of the variables (\r]i\ 2 }, which are the 
diagonal elements of the so-called x-matrix in the Oi 
basis (c/ the a, in (s], or the Xm.m in [H|). When 
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the error E acts on a short time we can expect the 
coefficients to be small and eq. © can be taken instead 
as a first order Taylor expansion of E, where the r]i 
with I > 1 play the role of a generator of the error. 
This idea was originally developed in [7J under the same 
setting, but aiming to study the generators directly led 
to limitations in the error model. Nevertheless, drawing 
the connection between the two opens the possibility 
of identifying a generator, which is essentially a Hamil- 
tonian with varying parameters r\i whose dynamics we 
can only observe on average through the {\m\ 2 )- This 
interpretation allows a different insight into the dynam- 
ics of the system, as many physical processes are better 
described by the action of a stochastic Hamiltonian 
rather than by an operator-sum representation arising 
from a system+bath picture. 

In conclusion, we have presented a method to char- 
acterize the spatial correlations occurring in a gate or 
process under study, showing its potential through a 
liquid state NMR QIP experiment. Second, we have 
pointed out the need of experimental feedback in order 
to arrive to a not only scalable but feasible protocol, as 
the one we introduced. Finally, we have analyzed the 
relevance of implementation errors, showing the need for 
strategies that are not only scalable but also robust. 

We thank J. Emerson, M. Silva and L. Viola for use- 
ful discussions, and J.S. Hodges and T. Borneman for 
invaluable assistance with the experiments. This work 
was supported in part by the National Security Agency 
(NSA) under Army Research Office (ARO) contract num- 
ber W911NF-05-1-0469. 



VI. APPENDIX A: ON THE CLIFFORD TWIRL 

Consider the twirling of the map § for m qubits in the 
Hilbert space Hm, of dimension 2 m . As shown in [l5j], a 
Haar twirl is equivalent to Clifford twirl as follows 

Pi = [ dUU^S(Up U^)U (10) 

JU(M) 

= wE4§(wi)c fc (ii) 

k=l 

where dU denotes the Haar measure on the unitary 
group U(2 m ). The sum on eq. (|11| runs over all the 
elements of the Clifford group for U(2 m ). In particular, 
this equation holds for a twirl on £/(2)® m , where the 
Ck are m-fold tensor products of the K — 24 elements 
of the Clifford group for one qubit, for the m qubits in 
Km- 



We can write the Clifford operators for one qubit 
as products SP of an clement S of the Symplec- 
tic group {exp(— ii/(ir/S)(a x + a y + er z )/\/3), v — 
0,1,2; exp(— ?(7r/4)(Tp), p — x,y,z} and an clement P 
of the Pauli group {I , a Xl a y , a z }. Due to some redun- 
dancy in the symplectic twirling, the minimum number of 
elements to achieve the 2-design property (eq. (0]) , as de- 
fined in [lj|) is actually half of the whole Clifford group: 
S can be taken from either Si = {exp(— iv(n /i){a x + 
<r v +(T z )/V3), i> = 0, 1,2} or S 2 = {exp(-i(7r/4)cr p ), p = 
x,y, z}. 

This number can be further decreased to K = 6 if we 
take an initial state po as described in the main text 
(eq. (3)) and we remark that we will be only interested 
in the projection of Pi^ onto given initial states of the 
form Pq = (I + a z )/2. There is then an additional 
redundancy in the values of these projections which 
enables us to chose P in the set {P\,P2\ with P\ = 1 

or er 2 and P 2 = cr x or a y . Here we have chosen p^ 
along z for definiteness (which is what we use in the 
experiment), but an equivalent result can be obtained 
for the x or y directions. Notice that the number 
of elements K required in eq. (|11[) depends on what 
quantity we will measure. In U(2)® m (as it is our case), 
whether some elements are redundant can be proved 
simply by straightforward calculation. 

Implementing a Haar twirl as defined on eq. (|10[) would 
require sampling over an infinite pool of random rotations 
Uk- The equivalence with the Clifford twirl allows us to 
implement the twirling with a finite set of gates: the 
K m Clifford operators in U(D). Nevertheless, although 
not infinite, the pool is of exponential size in m, thus 
again calling for a sampling strategy to implement the 
twirl approximately. Of course, for m small enough, it 
may be possible to implement the twirl exactly (as in our 
experiment, where m = 2 and the Clifford pool has size 
36). But on general grounds we can take on any of the 
two approaches. 

Depending on the experimental setup in question, one 
twirl may be more robust than the other (depending on 
what type of errors are expected for one-qubit gates), and 
also one pool may be easier to construct than the other. 



VII. APPENDIX B: COMBINING 7 (a) , AND 
7 (a,b) TQ RETRIEVE (\v2°b\ 2 ) 

To illustrate the mechanism of combining different 7's, 
consider measuring one and two qubits as follows. With 
Aj = {4Pj - 2)/3: 
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7 (a) = A((ICf>+E<IC/| 2 >+ E (K l j,k\) 2 + ---) and similarly 7 W. 



AaV b ({K°f)+ E(icji 2 >+ E (Ki, fe i 2 > + ---)+ 7vu(<k cc T> + E(iCi 2 )+ E <i<^i 2 > + 

^ j^a.b k>j^a,b ' ^ j^a.b k>j^a,b 



where . . . denote the corresponding higher order multi- 
body terms. If we prepare the qubits in a pure state, so 
V a = Vb = I, the combination 

7 W + 7 W _ 7 C<^) = 4 ^ 7 ^ |2) + g . |2) + j 

leaves only the collective coefficients involving both the 
qubits a and b. 

VIII. APPENDIX C: DETAILS ON THE 
EXPERIMENTAL IMPLEMENTATION 

The internal Hamiltonian of the system in the rotating 
frame is given by 

Hint - hJ2 ^*<P + » E ^f* { M k) (12) 

3=1 k>j=l 

where the chemical shifts, at our 9.4 T spectrometer, are 
of the order of kHz: ujs,i — 6650.6 Hz, = 1695.8 Hz, 
uj s ,3 = 4210.0 Hz, uj S j = -8796.7 Hz. The J-couplings 
are J 12 = 72.6 Hz, J 23 = 69.8 Hz and J u = 7.1 Hz, 
while J24 = 1.6 Hz, J13 = 1.3 Hz, and J34 = 41.6 Hz 
(according to the characterization of the sample we used; 
see also (l7|). 

The experimental initial state preparation over the 
4 qubits reported a correlation with the theoretical 
one that was on average 0.99 (0.98 the lowest). The 
correlation for the targeted qubits (a pair of qubits) was, 
in each case, between 1.00-0.99 (the pseudo-pure state 
preparation was designed to optimize this correlation). 

We implemented the twirl of pairs of qubits exactly 
using 36 Clifford operators. We had to pick one of the 8 
available 6-element subsets of Clifford operators for each 
qubit, and we chose the subset that performed best ex- 
perimentally: we applied each candidate to the thermal 
equilibrium state and compared the experimental perfor- 
mance with the theoretical one. This criterion coincided 
with choosing the subset that best took the equilibrium 
state to the 1/2 state for the qubit being twirled. 

To perform quantum state tomography (QST) on 4 
qubits with a liquid state NMR QIP, we used a set of 



18 readout pulses. Therefore the number of experiments 
required to measure one collective coefficient (l^ ^ 2 ) 
for a given pair of targeted qubits (a, b) for a particular 
gate under study was 648, plus 18 experiments to 
characterize the initial state |00)(//2)® 2 corresponding 
to preparing that pair in a pseudo-pure state. We 
performed QST of the full system therefore having a 
broader knowledge of the experimental performance, 
but this is not required by the protocol: only the target 
qubits must be measured. 

The negligibility of higher order multi-body terms 
in the gates under study is to be expected in liquid 
NMR QIP. In a simple model, these gates consist basi- 
cally of periods of free evolution of length r (the cor- 
responding propagator is U T = exp (— iHi nt T /K)) sep- 
arated by 7r-pulses on some of the qubits (so = 
exp (— i(7r/2 + e)cr£ ), already accounting for some error 
e). For example, the sequence for the gate Ie is: 

n3.4 x2 n3,4 sl,4 

t-it) x -T- it) x -t- ir) x - r - n) x - 

n3,4 n2 n 3,4 n1,4 

t - tt)_ x - t - ir)_ x -r- ir)^ x - t - n)_ x 

where r denotes free evolution for a time t, and tt)^ 
denotes a 7r-pulse (180 degree rotation) around the p-axis 
for the qubits j. 

Using the BCH formula [lj| is straighforward to see 
that in the building block VrlJ^, 3-body and higher or- 
der terms will appear with a factor at least -^j^ smaller 
respect to any possible 1-body and 2-body terms. For 
the values of Jj.kT involved in our experiment, we have 
-^2 < 0.14. This means that any possible 3-body and 
4-body terms would appear with a coefficient 10 times 
smaller than the ones for 1-body and 2-body terms. 

On the other hand, the simulation of the engineered 
pulse sequences used in the experiment showed that all 
the 3-body and 4-body terms appear with collective coef- 
ficients \Vj°k,f\ 2 , \Vj°lj>.k>\ 2 < °- 005 for thc CNOT and 
Ci2(0.4) gates, and < 0.002 for the rest. These are much 
smaller than the differences bctweend measured and pre- 
dicted values of |?7™i| 2 , which can be better explained 
as implementation errors in the protocol or genuine gate 
errors arising from imperfect refocusing. 
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